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HIGHLIGHTS 


•  We  show  a  particle-arrangement-preserving  method  for  sectioning  battery  electrodes. 

•  Energy-filtered  TEM  is  used  to  map  the  Li  distribution  in  half-charged  electrodes. 

•  Our  results  confirm  a  nucleation-limited  kinetic  model  for  LiFeP04  intercalation. 
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We  demonstrate  an  ex  situ  method  for  analyzing  the  chemistry  of  battery  electrode  particles  after 
electrochemical  cycling  using  the  transmission  electron  microscope  (TEM).  The  arrangement  of  particles 
during  our  analysis  is  the  same  as  when  the  particles  are  being  cycled.  We  start  by  sectioning  LiFeP04 
battery  electrodes  using  an  ultramicrotome.  We  then  show  that  mapping  of  the  Fe2+  and  Fe3+  oxidation 
state  using  energy-filtered  TEM  (EFTEM)  and  multivariate  statistical  analysis  techniques  can  be  used  to 
determine  the  spatial  distribution  of  Li  in  the  particles.  This  approach  is  validated  by  comparison  with 
scanning  transmission  X-ray  microscopy  (STXM)  analysis  of  the  same  samples  [Chueh  et  al.  Nanoletters, 
13  (3)  (2013)  866-72].  EFTEM  uses  a  parallel  electron  beam  and  reduces  the  electron-beam  dose  (and 
potential  beam-induced  damage)  to  the  sample  when  compared  to  alternate  techniques  that  use  a 
focused  probe  (e.g.  STEM-EELS).  Our  analysis  confirms  that  under  the  charging  conditions  of  the 
analyzed  battery,  mixed  phase  particles  are  rare  and  thus  Li  intercalation  is  limited  by  the  nucleation  of 
new  phases. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Modern  Li-ion  battery  systems  are  complex.  They  contain  com¬ 
posite  electrodes  with  a  size  distribution  of  active  particles  bound 
together  by  a  polymer  binder  in  a  complicated  3D  microstructure 
[1].  The  liquid  electrolyte  penetrates  between  the  particles  in  the 
electrode.  Particle  size,  porosity,  and  tortuosity  in  the  electrodes  are 
all  important  factors  to  consider  when  optimizing  the  performance 
and  degradation  of  a  battery  [2,3].  Ultimately,  understanding  how 
the  active  particles  are  arranged  in  an  electrode,  and  how  each  in¬ 
dividual  particle  undergoes  the  charge/discharge  reaction  will  yield 


*  Corresponding  author.  Tel.:  +1  925  294  1344;  fax:  +1  925  294  3231. 
E-mail  address:  jdsugar@sandia.gov  (J.D.  Sugar). 

0378-7753 /$  -  see  front  matter  ©  2013  Elsevier  B.V.  All  rights  reserved. 
http://dx.doi.Org/10.1016/j.jpowsour.2013.08.003 


a  fundamental  understanding  of  the  mechanisms  that  govern  the 
cycle  rates,  power  density,  energy  density,  and  efficiency  of  a  real 
cycling  battery.  Our  approach  is  to  disassemble  a  cycled  battery  and 
investigate  sections  of  the  post  mortem  electrode  particle-by¬ 
particle  to  understand  the  mechanisms  that  control  the  charge/ 
discharge  reactions.  In  previous  work,  we  have  relied  on  scanning 
transmission  X-ray  microscopy  (STXM)  to  chemically  analyze  the 
chemical  state  of  each  particle  in  a  LiFeP04  electrode  [4].  In  this 
study,  we  show  how  energy-filtered  TEM  data  can  be  used  to  probe 
chemical  state  in  the  same  system  at  higher  spatial  resolution.  The 
overall  good  agreement  between  the  two  approaches  validates  the 
Li  distribution  inferred  from  both  techniques  and  reinforces  the 
conclusion  of  the  previous  work:  that  under  some  conditions  the 
charging  rate  of  a  battery  can  be  determined  by  the  rate  of  nucle¬ 
ation  of  new  phases  within  each  particle. 
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LiFeP04  (LFP)  is  a  well-studied  battery  electrode  material  that  is 
attractive  because  of  its  relatively  low  cost,  flat  voltage  curve  at  3.5  V 
vs.  Li/Li+,  and  ease  of  safe  handling  [5,6].  When  cycled,  LFP  undergoes 
a  first-order  phase  transformation  from  the  Li-rich  LFP  phase  to  the 
Li-poor  FePC>4  (FP)  phase.  The  two  phases  have  very  limited  solid 
solubility  at  room  temperature,  and  both  phases  have  the  olivine 
structure  with  the  orthorhombic  unit  cell  [7].  There  has  been  much 
work  in  the  literature  investigating  the  mechanism  responsible  for 
the  Li-intercalation  reaction  in  these  materials  and  the  nature  of  the 
interphase  interface  during  the  transformation  [8-22].  An  important 
question  to  ask  is  whether  the  rate  of  the  LFP/FP  phase  trans¬ 
formation  is  limited  by  the  motion  of  the  interphase  boundary  or  the 
nucleation  of  the  new  phase.  If  the  phase  transformation  is  limited  by 
the  rate  at  which  the  phase  boundary  moves,  then  it  is  important  to 
differentiate  cases  that  proceed  by  different  mechanisms,  e.g.  a 
“core— shell”  [14]  or  “domino-cascade”  [11]  model.  Flowever,  if  the 
transformation  in  each  particle  is  rapid,  and  a  “many-particle”  or 
“mosaic”  pathway  [9,12,15]  better  describes  the  phase  trans¬ 
formation,  then  it  is  far  more  important  to  understand  how  to  in¬ 
crease  the  nucleation  rate  of  new  phases  during  cycling.  Despite 
many  studies,  some  very  recent,  different  groups  using  different 
techniques  reach  contradictory  conclusions  as  to  what  controls  phase 
transformations  in  LFP.  We  next  describe  how  the  mechanisms  can  be 
differentiated  by  measuring  the  fraction  of  particles  that  are  trans¬ 
forming  at  a  given  time.  We  then  show  the  successful  use  of  energy- 
filtered  TEM  to  perform  such  a  measurement,  which  unambiguously 
supports  nucleation-limited  kinetics. 

A  number  of  different  micro-  and  macroscopic  characterization 
techniques  have  been  used  to  probe  the  question  of  how  the  inter¬ 
calation  reaction  proceeds.  These  include  X-ray  and  neutron  diffrac¬ 
tion  [8,23-28],  transmission  electron  microscopy  (TEM) 
[10,13,21,29,30],  electrochemical  methods  [31,32],  and  recently,  X- 
ray  absorption  spectroscopy  [33,34]  and  X-ray  microscopy  [4].  The 
high  spatial  resolution  of  the  TEM  has  made  it  an  attractive  choice  for 
looking  at  the  micro  and  nanostructure  of  battery  electrode  materials. 
In  fact,  microscopy  can  distinguish  between  phase-boundary- 
migration  limited  and  nucleation-limited  kinetics  without  observing 
nucleation  events  or  measuring  the  time  it  takes  for  an  individual 
particle  to  change  phase.  The  key  distinguishing  parameter  is  the 
fraction  of  particles  in  a  mixed  charge  state  at  any  particular  time.  To 
see  this,  consider  the  simplified  charging  history  of  an  individual 
particle  shown  in  Fig.  1.  Initially  the  particle  is  uncharged.  At  time  to, 
the  charged  phase  nucleates.  The  charged  phase  grows,  so  that  the 
current  to  the  particle  is  a  constant  jo  (which,  of  course  depends  on  the 
charging  potential).  The  total  charging  time  of  the  particle  is 
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where  v  is  the  volume  and  ps  is  the  charge  density  per  unit  volume 
of  the  particle.  If  the  number  of  charging  particles  Nc  is  comparable 
to  the  total  number  of  particles  N,  than  the  total  charging  rate  J  will 
be  NJo  ~  Njo,  i.e.,  the  total  charging  rate  at  a  particular  voltage  will 
only  depend  on  how  fast  the  individual  particles  charge  after 
nucleating  a  new  phase  (schematically  shown  in  Fig.  1(b)).  All 
particles  are  charging  at  the  same  rate.  On  the  other  hand,  if  the 
number  of  particles  charging  at  any  particular  time  is  small,  then 

J  =  Ncj0  =  (NurAt)j0  =  Nurps  v  (2) 

where  Nu  is  the  number  of  untransformed  particles  and  r  is  prob¬ 
ability  per  unit  time  that  the  new  phase  nucleates.  A  substitution  of 
At  from  Equation  (1 )  into  (2)  removes  the  dependence  on  At.  Thus, 
the  charging  rate,  j,  is  determined  by  the  nucleation  rate  r  and  in¬ 
dependent  of  how  fast  an  individual  particle  charges  (schematically 
shown  in  Fig.  1(c)).  The  applied  potential  needed  to  charge  the 
battery  at  some  rate  depends  only  on  how  the  nucleation  rate 
varies  with  the  potential:  the  potential  needs  to  be  high  enough  to 
nucleate  the  number  of  charging  particles  per  unit  time  required  to 
give  the  specified  charging  rate.  This  model  neglects  such  compli¬ 
cating  effects  as  voltage  gradients  within  the  electrode  or  distri¬ 
butions  of  particle  sizes.  Nevertheless,  determining  if  the  fraction  of 
particles  in  a  mixed  state  of  charge  is  small  allows  one  to  determine 
whether  nucleation  or  individual  particle  charging  rates  determine 
the  overall  battery  charging  rate. 

A  number  of  electron  microscopy-based  techniques  have  been 
used  to  image  the  Li  distribution  in  battery  electrode  materials.  These 
techniques  have  the  potential  to  locate  and  count  the  number  of 
mixed-phase  particles  present.  They  include  electron  diffraction 
[10,29],  electron  energy  loss  spectroscopy  [13,35-39],  high- 
resolution  TEM  [21],  and  annular  bright-filed  imaging  [30].  For 
many  of  these  techniques,  it  is  necessary  to  take  particles  and  disperse 
them  on  a  TEM  support  grid,  which  destroys  the  native  microstructure 
that  existed  in  an  operating  battery.  Our  goal  was  to  develop  a  tech¬ 
nique  that  would  allow  us  to  observe  active  LFP  particles  in  the  same 
arrangement  as  they  exist  in  a  battery  during  cycling,  and  to  image  the 
distribution  of  Li  in  each  particle  of  the  electrode  after  the  battery  was 
cycled  to  a  specific  state  of  charge  to  determine  the  fraction  of  mixed 
phase  particles  during  charging.  We  accomplished  this  using  a  com¬ 
bination  of  ultramicrotomy  for  sample  preparation  (see  for  example 
Ref.  [40]),  and  energy-filtered  TEM  [41  ].  Our  techniques  are  advanta¬ 
geous  for  two  reasons.  1 )  The  arrangement  of  particles  in  the  battery  is 
preserved  during  analysis  because  sections  are  microtomed  from  a 
real  electrode.  2)  Also,  we  use  a  parallel  electron  beam  in  energy- 
filtered  TEM  rather  than  a  focused  probe  as  in  STEM— EELS  (scan¬ 
ning  transmission  electron  microscopy-electron  energy  loss  spec¬ 
troscopy),  which  reduces  some  effects  of  electron  beam  damage. 

We  compared  our  results  to  scanning  transmission  X-ray  mi¬ 
croscopy  (STXM)  data  and  found  good  agreement,  giving  us  confi¬ 
dence  in  the  analysis.  Ultimately,  our  ex  situ  forensic  analysis  on  LFP 
electrodes  under  the  charging  conditions  of  this  study  suggests  that 
the  nucleation  of  new  phases  is  the  most  important  rate¬ 
determining  step  in  the  phase  transformation. 

2.  Experimental  procedures 

2  A.  Electrochemical  cycling 


Fig.  1.  (a)  Schematic  showing  the  charging  of  an  individual  LiFeP04  particle.  The  fully 
lithiated  particle  has  a  volume,  v,  lithium  concentration,  ps,  and  charges  in  time  At.  The 
transformation  rate  is  limited  by  phase  boundary  migration  and  nucleation  in  (b)  and  (c), 
respectively. 


Standard  coin-cell  geometry  batteries  were  prepared  with  a  Li 
anode  and  an  approximate  35-pm  thick  carbon-coated  LFP / 
graphite/PVDF  composite  cathode.  The  average  LFP  particle  size 
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was  ~  220  nm,  and  in  general,  each  particle  of  LFP  is  a  single  crystal. 
The  electrolyte  was  1.2  M  LiPF6/ethylene  carbonate  (EC)/ethyl 
methyl  carbonate  (EMC)  with  a  3:7  weight  ratio  of  EC: EMC.  The 
battery  was  electrochemically  cycled  five  times  to  obtain  a  stable 
capacity  of  150  mA  h  g-1.  We  then  charged  a  battery  at  a  rate  of 
1  C— 75  mA  h  g_1  (50%  state  of  charge).  Electrochemical  data  can  be 
found  in  the  supplemental  information  of  Ref.  [4].  After  reaching 
the  desired  state  of  charge,  the  cell  was  rapidly  disassembled  and 
rinsed  with  excess  dimethyl  carbonate  in  a  dry  room.  The  disas¬ 
sembly  procedure  was  completed  within  4  min  of  terminating  the 
charge,  which  is  faster  than  typical  voltage  relaxation  times  [42,43]. 
The  LFP  particles  are  mostly  not  in  contact  with  one  another  after 
the  electrolyte  is  removed  [1],  which  prevents  further  Li  redistri¬ 
bution  after  disassembly.  It  is  possible,  however,  for  Li  to  redis¬ 
tribute  within  a  single  particle,  however,  we  will  show  that  very 
few  of  the  particles  we  measured  exhibit  two-phase  coexistence 
within  a  single  particle. 

2.2.  Transmission  electron  microscopy 

After  the  battery  coin  cells  were  disassembled,  the  LFP  elec¬ 
trodes  were  removed.  Two  methods  were  used  to  examine  indi¬ 
vidual  particles  of  LFP  by  TEM  after  being  electrochemically  cycled 
in  the  coin  cell.  In  the  first  method,  LFP  particles  were  scraped 
from  the  electrode  with  a  razor  blade.  These  particles  were 
dispersed  in  ethanol  and  ultrasonicated  for  approximately  1  h.  The 
particles  were  pipetted  onto  a  Cu  TEM  support  grid  (Ted  Pella  Inc.). 
For  the  second  method,  a  small  piece  of  the  electrode  was  cross- 
sectioned  into  electron  transparent  slices  using  the  ultramicro¬ 
tome  (Leica  EM  UC7).  The  ultramicrotome  produced  thin  slices 
across  the  entire  width  of  the  electrode.  The  operation  of  the  ul¬ 
tramicrotome  to  produce  thin  cross  sections  is  shown  schemati¬ 
cally  in  Fig.  2(a).  Slices  of  electrode  are  made  in  which  the 
arrangement  of  active  LFP  particles  in  each  slice  is  not  disturbed 
relative  to  their  arrangement  in  the  coin  cell.  The  thin  cross- 
sections  spanned  the  entire  width  of  the  electrode  from  the 
electrode/electrolyte  interface  to  the  current  collector  interface.  An 
optical  micrograph  of  a  cross  section  prepared  with  the  ultrami¬ 
crotome  is  shown  in  Fig.  2(b).  The  entire  35-pm  thick  LFP  electrode 
is  visible,  and  a  portion  of  the  Al  current  collector  has  delaminated 
from  the  electrode  during  the  sectioning. 

Transmission  electron  microscopy  was  performed  using  a  JEOL 
2010F  operating  at  200  kV  equipped  with  a  Gatan  860  GIF  EELS 
spectrometer.  In  addition  to  imaging  the  particles,  EELS  (electron 
energy  loss  spectroscopy)  spectrum  imaging  was  performed  in 
which  an  EELS  spectrum  is  recorded  at  each  pixel  of  an  image.  Two 
acquisition  techniques  were  used  to  collect  EELS  spectrum  images 
and  these  are  depicted  in  Fig.  3.  In  Fig.  3(a),  a  focused  STEM 
(scanning  transmission  electron  microscopy)  probe  is  used  to 
collect  an  EELS  spectrum  at  each  pixel.  One  spectrum  is  collected, 
the  probe  is  moved,  and  then  the  next  spectrum  is  collected.  The 
other  technique,  EFTEM  (energy- filtered  TEM)-spectrum  imaging, 
is  shown  in  Fig.  3(b).  This  technique  utilizes  a  conventional,  parallel 
electron  beam  and  collects  an  image  at  a  particular  energy  loss  with 
an  energy-selecting  slit.  Then,  the  energy-selecting  slit  is  moved  an 
energy,  A E,  and  another  image  is  collected.  In  this  way,  an  EELS 
spectrum  is  constructed  at  each  pixel,  and  the  same  information 
can  be  collected  as  in  the  STEM-EELS  case.  We  collected  data  with 
an  energy  slit  width  of  5  eV  and  a  AE  of  1  eV  for  each  energy  plane. 
Usually,  there  are  more  data  artifacts  introduced  with  EFTEM  data 
collection  than  STEM-EELS.  These  artifacts  cannot  be  corrected 
during  acquisition  easily  (e.g.  dark  current  drift,  ghosting  on  the 
CCD,  spatial  drift,  energy  drift).  We  collected  our  data  with  the 
SmartEFTEM  acquisition  package  [44],  which  eliminated  many  is¬ 
sues  associated  the  CCD  during  data  collection. 
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Fig.  2.  (a)  Schematic  showing  the  operation  of  a  microtome  to  fabricate  cross  sections 
of  battery  electrodes.  The  cross  sections  are  electron  transparent,  and  they  span  the  full 
width  of  the  battery  electrode.  An  optical  micrograph  of  a  battery  cross  section  is 
shown  in  (b).  The  current  collector  has  delaminated  from  the  electrode  during 
sectioning.  The  electrode  spans  across  holes  in  the  TEM  support  grid,  where  the  TEM 
can  resolve  each  individual  particle. 

2.3.  Data  processing  and  analysis 

Several  artifacts  are  introduced  during  EFTEM-SI  data  acquisi¬ 
tion  that  must  be  corrected  for  such  as  non-isochromaticity,  spatial 
drift,  and  energy  drift.  Correction  of  these  artifacts  post-acquisition 
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Fig.  3.  A  schematic  illustrating  the  difference  between  STEM-EELS  (a)  and  EFTEM  (b) 
spectrum  imaging.  For  STEM-EELS,  a  focused  STEM  probe  is  scanned  across  a  region  of 
the  sample  and  EELS  spectra  are  collected  at  each  pixel.  EFTEM  spectrum  image  data  is 
collected  with  a  parallel  electron  beam,  and  images  are  acquired  of  each  energy  plane. 
The  end  result  in  both  cases  is  data  at  each  ( x ,  y,  energy)  triplet. 


J.D.  Sugar  et  al.  /  Journal  of  Power  Sources  246  (2014)  512-521 


515 


is  described  in  detail  in  the  literature  [45].  We  corrected  for  non- 
isochromaticity  using  the  position  of  the  zero-loss  peak,  and  we 
used  the  freely  available  SDSD  correction  for  spatial  drift  [46].  We 
did  not  correct  for  energy  drift  because  we  determined  that  the 
total  energy  drift  during  acquisition  was  much  smaller  than  the 
width  of  the  energy-selecting  slit  (5  eV),  and  therefore  it  would  not 
make  a  significant  difference  in  the  data  analysis. 

To  determine  the  Li  distribution  from  EFTEM-SI  data,  we 
examined  the  state  of  charge  of  the  Fe  from  the  EELS  spectra.  The 
literature  has  shown  that  the  position  of  the  Fe  L  edge  shifts  by 
approximately  2  eV  when  Fe  changes  from  Fe2+  in  LiFeP04  to  Fe3+ 
in  FeP04  [13].  This  shift  is  difficult  to  measure  with  EFTEM-SI 
because  the  width  of  the  energy-selecting  slit  is  5  eV  and  a 
typical  energy  resolution  of  the  EELS  spectrometer  (FWHM  of  the 
zero-loss  peak)  is  1.5  eV.  While  it  is  still  possible  to  measure  the  Fe  L 
edge  shift,  this  task  was  made  easier  by  our  use  of  multivariate 
statistical  analysis  (MSA)  techniques.  MSA  describes  a  series  of 
mathematical  techniques  used  for  matrix  decomposition  and  is 
commonly  used  in  analytical  electron  microscopy  [47-49].  A 
spectrum  image  data  cube,  D,  of  p  channels  and  dimensions  (x-y) 
can  be  represented  by  the  more  compact  form: 

D=CSr  (3) 

where  C  is  a  matrix  of  spatial  pixel  weightings  of  dimension  (x  y )  •  n, 
and  S  represents  the  p  ■  n  matrix  of  “pure”  component  spectra.  This 
technique  assumes  that  the  data  follows  a  linear  additive  model.  C 
and  S  can  be  calculated  from  a  variety  of  MSA  techniques  including 
principal  component  analysis  (PCA),  singular  value  decomposition 
(SVD),  and  multivariate  curve  resolution  (MCR).  MCR  has  been 
shown  to  be  successful  in  the  analysis  X-ray  absorption  spectra 
from  battery  electrode  materials  [50]. 

We  used  a  “spatially  simple”  multivariate  curve  resolution  al¬ 
gorithm  developed  at  Sandia  (AXSIA)  that  has  been  described  in 
detail  in  the  literature  [47,48,51-55].  The  results  of  the  MSA  anal¬ 
ysis  yielded  pure  spectral  components  that  corresponded  to  the 
Fe2+  and  the  Fe3+  states  in  addition  to  other  carbon  and  spectral 
background  components.  Reference  spectra  for  the  pure  Fe2+  and 
Fe3+  states  could  be  extracted  from  the  raw  data  by  locating  the 
most  intense  pixels  of  the  pure  components  in  the  MSA  analysis. 
These  reference  spectra  were  also  used  for  multiple  linear  least 
squares  (MLLS)  fitting  to  map  the  Fe  state  of  charge  in  an  area  of  the 
sample.  MLLS  also  solves  Equation  (3),  but  instead  of  calculating  S 
from  MSA,  the  spectra  used  in  S  are  specified  by  the  user  and  are 
extracted  directly  from  the  data  or  selected  from  a  library  of 
reference  spectra.  Then,  a  simple  matrix  division  solves  for  the 
coefficients  in  C.  The  Fe  state-of-charge  map  was  determined  by: 

Fe  state  of  charge  =  - - ^FeP°4 -  (4) 

*LiFeP04  +  'FeP04 

where  the  intensities,  /uFeP04  and  fFeP04  are  determined  by  the  co¬ 
efficients  in  C  of  the  MLLS  fit.  When  the  Fe  state-of-charge  is 
calculated  to  be  0  at  a  particular  pixel,  the  Fe  is  in  the  Fe2+  state  and 
Li  is  present,  when  it  is  1,  it  is  Fe3+.  The  determination  of  the  Li 
distribution  after  cycling  of  the  LFP  electrode  material  allowed  us  to 
learn  about  the  mechanisms  of  the  phase  transformation  during 
electrochemical  cycling. 

3.  Results 

The  first  step  in  determining  the  Li  distribution  in  cycled  battery 
electrode  material  was  to  determine  the  spectral  reference  shapes 
for  Fe2+  in  LiFeP04  and  Fe3+  in  FeP04.  It  is  well  known  from  the 
literature  that  our  monochromated  system  can  easily  differentiate 


between  the  two  spectral  shapes  [13].  It  was  not  clear  that  a  non- 
monochromated  system  would  be  able  to  differentiate  so  easily 
because  the  energy  resolution  ( ~  1.5  eV)  is  roughly  the  same  as  the 
Fe  L  edge  shift.  Reference  shapes  for  FeP04  was  determined  in  two 
ways:  1)  The  chemical  delithiation  of  LiFeP04  by  soaking  in  a  so¬ 
lution  of  N02BF4  in  acetonitrile,  and  2)  the  electrochemical  deli¬ 
thiation  of  LiFeP04  using  a  coin  cell  as  described  in  the 
Experimental  procedures  above.  The  resulting  EELS  spectra  for  both 
experiments  is  shown  in  Fig.  4(a)  and  (b). 

Regardless  of  whether  the  FeP04  is  prepared  by  chemical  deli¬ 
thiation  or  electrochemical  delithiation,  the  spectra  show  a  shift 
1.5—2  eV  for  the  peak  position  of  the  Fe  L3  edge.  This  peak  shift  is 
observable  in  our  non-monochromated  system,  which  is  somewhat 
surprising  because  the  energy  resolution  measured  by  the  FWHM 
of  the  zero-loss  peak  is  typically  approximately  1.5  eV.  However, 
these  spectra  show  that  it  is  possible  to  use  STEM-EELS  to  map  the 
Fe  state-of-charge.  Unfortunately,  the  use  of  a  focused  electron 
probe  in  STEM  mode  causes  beam  damage  that  ultimately  limits 
the  usability  of  STEM. 

The  effects  of  beam  damage  are  shown  in  Fig.  4(c).  Here,  a  STEM- 
probe  was  kept  stationary  on  a  single  FeP04  particle  and  an  EELS 
spectrum  was  recorded  every  second  (with  a  1  s  exposure)  for  50  s. 
There  are  two  visible  artifacts  that  result  from  beam  damage.  The 
HAADF  image  shows  a  small  hole  in  the  particle  at  the  location  of 
the  focused  electron  beam  (marked  with  an  arrow).  So,  the  focused 
electron  beam  has  removed  material  and  caused  a  large  decrease  in 
the  available  signal  for  EELS  as  shown  by  the  difference  in  the  EELS 
intensities.  Also,  there  is  a  visible  shift  in  the  position  of  the  Fe  L3 
edge  of  1  eV,  where  the  last  spectrum  recorded  has  a  peak  position 
1  eV  lower  than  the  first  spectrum.  This  is  consistent  with  the 
beam-induced  reduction  of  Fe  from  Fe3+  to  Fe2+  and  evidence  in 
the  literature  that  the  electron  beam  can  cause  the  reduction  of  Fe 
ions  [56].  Electron-beam  induced  reduction  of  cations  is  a  common 
occurrence  in  the  electron  microscope  [57,58]. 

While  it  is  unlikely  that  one  would  need  a  50  s  exposure  to 
measure  Fe  EELS  spectra  during  a  spectrum  image,  this  time  is  used 
to  illustrate  the  extreme  case,  the  onset  of  beam-induced  Fe 
reduction  can  occur  at  much  shorter  times.  The  results  of  Pan  et  al. 
[56]  suggest  that  lowering  the  accelerating  voltage  of  the  micro¬ 
scope  to  100  kV  has  the  potential  to  slow  the  damage,  however,  it  is 
unclear  whether  the  damage  is  caused  by  direct  atomic  displace¬ 
ment  or  secondary  events.  Even  if  the  primary  electron  beam  en¬ 
ergy  is  reduced,  it  may  be  necessary  to  increase  the  acquisition  time 
for  high  quality  spectra.  Since  the  total  electron  fluence  (electrons/ 
nm-2)  determines  the  onset  of  electron-beam  induced  Fe  reduc¬ 
tion,  lowering  the  beam  energy  might  not  prevent  sample  damage. 
Therefore,  the  use  of  STEM-EELS  for  mapping  the  Fe  state-of- 
charge  in  LiFeP04-based  materials  will  be  limited,  particularly  in 
our  characterization  system  because  the  measurement  itself  could 
change  the  Fe  state  of  charge.  Therefore,  we  determined  that 
EFTEM-SI  would  be  a  more  successful  approach  because  one  can 
spread  the  beam  to  a  larger  area,  greatly  decreasing  the  electron 
dose  to  the  area  of  interest  and  minimizing  the  effects  of  beam- 
induced  artifacts. 

The  use  of  EFTEM-SI  and  MSA  to  determine  the  Fe  state  of 
charge  is  shown  in  Fig.  5.  A  zero-loss  image  of  LiFeP04  particles 
electrochemically  charged  to  50%  state  of  charge  and  distributed  on 
a  Si3N4  support  grid  is  shown  in  Fig.  5(a).  The  scree  plot  used  to 
generate  the  number  of  principal  components  is  shown  in  Fig.  5(b). 
The  first  5  components  deviate  significantly  from  the  linear  fit  to 
the  random  noise  components  (component  #6  and  above).  The 
“score”  (pixel  weighting)  images  for  component  #2  and  component 
#4  show  that  these  correspond  to  the  Si3N4  support  that  holds  the 
particles.  The  spectra  and  score  images  for  components  #  1, 3,  and  5 
are  shown  in  Fig.  5(c).  The  spectra  for  component  #3  and  #5  show  a 
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Fig.  4.  EELS  spectra  from  reference  materials  of  LiFeP04  and  FeP04  prepared  by  (a)  chemical  delithiation  of  LiFeP04  in  N02BF4,  and  (b)  electrochemical  delithiation  of  LiFeP04  in  a 
battery  coin  cell.  In  (c),  the  results  of  collecting  an  EELS  spectrum  every  second  for  50  s  is  shown.  There  is  a  reduction  in  the  Fe  L  edge  intensity  and  a  shift  of  the  peak  position 
consistent  with  the  reduction  of  Fe.  The  HAADF  STEM  image  shows  that  the  focused  electron  beam  has  removed  material  and  made  a  hole  (marked  with  arrow). 


shift  in  the  Fe  L3  edge  of  approximately  3  eV.  Although  the  detailed 
fine  structure  differences  (e.g.  Ref.  [13])  are  not  visible  with  an 
energy  step  of  1  eV  in  the  EFTEM  spectra,  the  observed  peak  shift  is 
consistent  with  the  observed  shift  in  position  of  the  Fe  L3  edge 
when  Fe  changes  from  a  3+  to  2+  oxidation  state.  Therefore, 
components  #3  and  #5  are  consistent  with  LiFePCH  and  FePCH, 
respectively.  Component  #1  causes  most  of  the  centers  of  the 
particles  to  light  up,  and  since  most  of  these  particles  are  round,  we 
expect  them  to  be  thicker  in  the  middle.  The  thickness  map  is  also 
consistent  with  component  #1  identifying  regions  where  the 
relative  thickness  (t/A)  is  greater  than  approximately  1.5.  Therefore, 
component  #1  is  identifying  sample  region  that  is  thick,  which 
causes  non-linear  effects  in  the  spectra  and  a  peak  position  of  the 
L3  edge  that  is  mixed  between  the  apparent  position  for  2+  and  3+ 
Fe.  When  samples  are  thick,  plural  scattering  occurs  and  causes  a 
redistribution  of  intensity  to  higher  energy  losses  [59-61].  This 
effect  can  cause  the  apparent  mixing  of  Fe2+  and  Fe3+  states  in  the 
fitting  because  intensity  is  artificially  higher  at  energies  between 
the  L3  edge  peak  intensities  for  the  Fe2+  and  Fe3+  reference  states. 

The  EFTEM-SI  measurements  of  Fe  oxidation  state  were  vali¬ 
dated  by  a  comparison  with  measurements  obtained  using  scan¬ 
ning  transmission  X-ray  microscopy  (STXM).  Ultimately,  it  was 
desired  to  validate  the  results  of  the  EFTEM-SI  data  with  another 
measurement  of  Fe  oxidation  state.  STXM  spectra  are  similar  to 
EELS  spectra  in  that  they  both  measure  absorption  of  the  incoming 
illumination  by  the  sample.  However,  STXM  uses  illumination 
generated  from  a  synchrotron  X-ray  source  (Advanced  Light  Source 


at  LBNL),  while  EELS  uses  electron  illumination  in  the  TEM.  The 
advantages  of  the  STXM  are  the  increased  energy  resolution  and 
decreased  sensitivity  of  the  sample  to  X-ray  beam  damage,  but  the 
spatial  resolution  of  STXM  is  only  ~30  nm,  much  less  than  what  is 
possible  in  the  TEM.  Fig.  6  shows  the  same  particles  supported  on  a 
Si3N4  grid  measured  with  both  STXM  and  TEM.  The  details  of  the 
STXM  can  be  found  in  Ref.  [4].  It  is  much  easier  to  identify  the  in¬ 
dividual  particle  boundaries  and  when  particles  overlap  in  the  TEM 
zero-loss  image,  illustrating  the  higher  spatial  resolution  of  the 
TEM. 

The  Fe  state-of-charge  map  in  Fig.  6  is  shown  in  two  ways  that 
can  be  compared  to  the  STXM  data.  Both  the  output  of  the  MSA 
analysis  (similar  to  Fig.  5)  and  the  Fe  state  of  charge  calculated  by 
MLLS  (Equation  (4))  is  shown  and  both  are  relatively  consistent 
with  the  STXM  state-of-charge  map.  The  MSA  output  shows  regions 
that  are  mostly  red  and  mostly  green,  and  these  regions  are 
consistent  with  the  red  and  green  locations  in  the  STXM  data. 
However,  the  particle  centers  are  not  identified  as  either  red  or 
green  in  the  MSA  output  because  plural  scattering  and  non-linear 
thickness  effects  create  another  component  for  thicker  regions  (as 
seen  in  Fig.  5).  Also,  the  MLLS  calculation  results  in  many  of  the 
particle  centers  being  yellow  (a  linear  combination  of  the  red  2+ 
and  green  3+  states).  This  is  also  caused  by  plural  scattering  and 
non-linear  thickness  effects  in  the  EELS  spectra.  In  order  to  obtain 
better  Fe  state-of-charge  maps  with  fewer  artifacts  from  the  non¬ 
linear  thickness  effects,  we  prepared  samples  with  the  ultrami¬ 
crotome.  This  results  in  samples  with  a  more  uniform  thickness 
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Fig.  5.  (a)  Zero-loss  image  of  50%  state  of  charge  particles  on  a  SisN4  support  grid,  (b)  Scree  plot  generated  by  spatial  simplicity  MCR  decomposition  of  the  EFTEM  data  set.  The  pure 
component  spectra  (loadings)  and  their  corresponding  score  images  are  shown  in  (c)  for  components  1,3,  and  5.  These  components  correspond  to  particles  and  are  consistent  with 
LiFeP04  in  component  5,  FeP04  in  component  3,  and  thick  material  in  which  non-linear  thickness  effects  dominate  in  component  1. 


across  the  sample  and  fewer  thick  regions  that  have  harder  to 
interpret  EELS  spectra. 

Fig.  7  shows  that  when  EFTEM-SI  analysis  is  performed  on 
ultramicrotomed  samples,  there  is  good  agreement  with  the  STXM 
data  and  fewer  thickness  artifacts  are  observed.  A  comparison  is 
shown  for  three  sample  regions.  In  each  region,  the  Fe  state  of 
charge  is  measured  with  both  STXM  and  EFTEM-SI.  The  zero-loss 


TEM  images  in  Fig.  7(a)  show  that  the  spatial  resolution  of  the 
TEM  is  sufficient  to  determine  when  particles  overlap.  Fig.  7(b) 
shows  the  STXM  maps  of  the  Fe  oxidation  state.  In  the  STXM  maps, 
regions  yellow  in  color  suggest  that  the  Fe  has  a  mixed  oxidation 
state.  It  is  not  possible  to  tell  from  the  STXM  data  alone  whether  the 
mixed  oxidation  state  occurs  because  an  individual  particle  has 
both  the  LiFePCU  and  FePCU  phases,  or  because  two  particles  with 
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Fig.  6.  Comparison  of  EFTEM  and  STXM  data  for  LiFeP04  particles  suspended  on  a  SbN4  support  film.  A  comparison  of  the  TEM  MSA  results  and  the  STXM  results  show  that  the  red 
regions  and  green  regions  are  identical.  The  non-linear  thickness  effects  cause  the  MLLS  analysis  on  the  EFTEM  data  to  have  more  yellow  than  the  STXM  data.  (For  interpretation  of 
the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


different  oxidation  states  are  overlapping  each  other.  The  TEM  data 
allows  us  to  more  easily  discern  when  particle  overlap  is  occurring. 

The  results  of  the  TEM  analysis  suggest  that  the  particles  are 
generally  either  in  a  Fe2+  or  a  Fe3+  state  with  mixed  oxidation  states 
only  occurring  when  there  is  significant  particle  overlap.  The  results 
of  MSA  calculation  are  shown  in  Fig.  7(c).  Here,  only  the  components 
for  the  pure  Fe2+  (red)  and  Fe3+  (green)  are  shown,  the  components 
for  the  C  support  and  background  have  been  omitted.  There  are  still 
some  particles  that  are  black  in  the  center  because  of  thickness  ef¬ 
fects,  but  the  MSA  allows  us  to  clearly  identify  regions  that  exhibit 
the  two  different  oxidations  states  as  shown  by  the  - 1.5  eV  shift  of 
the  Fe  L3  edge  in  the  component  spectra.  The  MSA  data  is  used  to 
identify  the  most  intense  pixels  for  the  Fe2+  and  Fe3+  states.  Then, 
spectra  from  these  regions  are  used  as  reference  shapes  in  the  MLLS 
fit  to  create  the  Fe  state  map  (Equation  (4))  in  Fig.  7(d).  A  comparison 
of  the  STXM  maps  in  Fig.  7(b)  and  the  Fe  maps  determined  from 
MLLS  in  Fig.  7(d)  are  generally  in  good  agreement.  That  two  such 
different  analysis  techniques  give  the  same  result  gives  considerable 
confidence  that  the  charge  maps  are  correct. 

4.  Discussion 

As  described  in  the  introduction,  to  answer  the  question  of 
whether  charging  is  determined  by  nucleation  or  individual  par¬ 
ticle  charging  (phase  boundary  migration)  rates  one  needs  to 
determine  the  fraction  of  mixed  phase  particles  at  any  time.  The 
main  experimental  difficulty  in  doing  this  is  to  determine  whether 
the  “yellow”  regions  in  Fig.  7  correspond  to  phase  boundaries  in 
individual  particles  or  whether  they  are  created  by  overlap  be¬ 
tween  two  independent  particles.  Consistent  with  the  STXM 
chemical  analysis,  the  results  of  the  EFTEM  analysis  suggest  that 
particles  are  rarely  in  a  mixed  state  of  charge.  Indeed,  the  high- 
spatial  resolution  of  the  TEM  allows  us  to  definitively  determine 
that  many  particles  with  a  mixed  state  of  charge  appearance  are 
actually  overlapping  other  particles  in  the  projected  view.  We  now 
discuss  some  of  the  subtleties  associated  with  this  determination 


and  how  high-resolution  data  and  the  combination  of  EFTEM  and 
STXM  can  be  useful  in  resolving  them. 

Three  regions  in  Fig.  7(d)  have  been  marked  (a,  b,  and  c)  for 
further  investigation.  In  a  it  is  ambiguous  whether  there  are  two 
overlapping  particles  or  a  single  particle  that  exhibits  both  Fe  states. 
A  higher-magnification  view  of  this  region  is  shown  in  Fig.  8(a).  The 
particle  edges  have  been  outlined  for  clarity  with  a  white  dotted 
line.  This  image  shows  that  overlapping  particles  causes  the 
appearance  of  mixed  Fe  states  in  the  projection  view.  Similarly,  in 
Fig.  8(b),  there  are  several  overlapping  particles,  and  the  STXM  data 
does  not  show  one  of  the  green  particles.  The  spatial  resolution  of 
the  STXM  makes  it  difficult  to  discern  particle  boundaries,  and  this 
particle  could  be  too  thin  to  give  a  sufficient  STXM  signal.  It  is 
possible  that  the  red  region  seen  in  STXM  is  entirely  due  to  the 
overlapping  red  particle  seen  in  EFTEM.  Thus  there  is  no  conclusive 
evidence  for  a  mixed  phase  particle  in  this  region. 

The  situation  in  Fig.  8(c)  is  more  complicated.  In  EFTEM,  there 
are  two  central  particles  that  show  a  significant  amount  of  mixed 
phase  character.  The  particle  in  the  lower  left  is  red  on  the  right  and 
transitions  to  green  on  the  left.  The  BF  TEM  image  reveals  that  a 
small  particle,  circled  in  the  image,  overlaps  the  upper  right  portion 
of  the  large  particle.  This  overlap  is  likely  the  cause  of  the  mixed 
phase  character  in  this  lower  particle.  The  large  upper  particle 
overlaps  other  particles  and  contains  a  fairly  high  mixed  charge 
component  (yellow)  distributed  uniformly  throughout  the  particle 
in  EFTEM,  while  STXM  suggests  a  phase  boundary  on  the  upper 
right  of  the  particle.  A  measurement  of  the  relative  thickness  of  this 
particle  with  EELS  shows  that  t/A  is  as  high  as  2.5.  Since  this  particle 
is  particularly  thick,  it  is  possible  that  non-linear  thickness  effects 
on  the  EELS  spectra  due  to  plural  scattering  cause  the  spurious 
appearance  of  mixed  states  of  charge  as  discussed  above.  Another 
possibility  is  that  the  particle  has  rotated  slightly  between  the 
STXM  and  EFTEM  measurements.  Then,  the  “red”  charged  region 
detected  by  STXM  could  be  underneath  the  “green”  region  causing 
“yellow”  in  the  projected  view  of  the  EFTEM  state  of  charge  map. 
Certainly,  STXM  is  less  sensitive  to  thickness  effects  than  TEM  and  is 
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Fig.  7.  Comparison  of  EFTEM  and  STXM  Fe  state-of-charge  mapping.  Bright-field  images  of  electrode  particles  are  shown  in  (a).  The  Fe  oxidation  state  measured  from  STXM  is 
shown  in  (b).  This  is  compared  to  the  Fe  state  measured  by  EFTEM  in  (c)  and  (d).  In  (c),  the  LFP  and  FP  components  of  the  MSA  analysis  are  mapped.  Good  agreement  is  found 
when  compared  to  the  STXM  data.  In  (d),  intense  LFP  and  FP  regions  are  identified  from  (c),  and  reference  spectra  shapes  are  extracted  and  used  to  perform  the  MLLS  fit  in  (d). 
Typical  spectral  shapes  for  the  two  techniques  are  shown  in  (e).  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this 
article.) 


less  likely  to  have  the  artificial  appearance  of  mixed  phase  particles, 
but  STXM  cannot  discern  the  particle  boundaries  easily  at  di¬ 
ameters  of  ~220  nm.  Despite  these  subtle  differences  in  the  state 
of  charge  maps,  the  general  conclusion  that  a  small  fraction  of  total 
particles  is  actively  charging  at  a  given  time  is  supported  by  the 
results  of  both  techniques.  It  is  clear  that  a  combination  of  the  two 
techniques  provides  more  information  and  confidence  for  mapping 
charge  states  and  is  better  than  one  measurement  alone. 

Analysis  of  a  large  number  of  particles  in  the  STXM  for  which 
there  was  little  ambiguity  about  particle  overlap  from  the  same 
sample  shows  that  particles  exhibit  both  Fe  oxidation  states  only 
approximately  2%  of  the  time  [4].  The  other  98%  are  either 
completely  Fe2+  or  Fe3+.  Mixed  phase  particles  are  certainly  rare.  A 


recent  study  using  in  situ  X-ray  measurements  found  that  during 
the  intercalation  reaction  LFP  and  FP  particles  exist  with  a  Li 
composition  within  the  miscibility  gap  [22,62].  This  measurement 
confirms  that  nucleation  is  a  critical  rate-limiting  step  in  the  re¬ 
action  because  the  existence  of  metastable  states  implies  the  dif¬ 
ficulty  of  nucleation.  Such  a  measurement  on  ensembles  of 
particles  sheds  insight  into  the  nucleation  process  itself,  however, 
to  be  more  quantitative  about  exactly  how  many  particles  are  in  a 
mixed  state  of  charge  requires  a  particle-by-particle  measurement 
with  careful  analysis  of  particle  boundaries  at  high  spatial  resolu¬ 
tion.  Then,  particle  overlap  can  be  identified.  This  is  not  possible  to 
infer  from  STXM  data  alone,  and  requires  a  high-resolution  tech¬ 
nique  like  we  show  here  with  the  TEM.  The  motivation  for  such  an 
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C)  500  nm 


Fig.  8.  High-magnification  views  of  regions  marked  a,  b,  and  c  in  Fig.  7(d)  with  corresponding  bright-field  (BF)  TEM  images.  Some  particle  boundaries  are  outlined  in  the  BF  TEM 
images  for  clarity.  Close  inspection  shows  that  for  these  particles,  mixing  of  Fe2+  and  Fe3+  can  be  explained  as  a  result  of  overlapping  particles  in  projection.  In  general,  we  do  not 
observe  particles  exhibiting  multiple  states  of  charge.  Some  of  the  subtleties  of  interpretation  are  discussed  in  the  text. 


analysis  is  that  the  fraction  of  particles  that  are  in  the  process  of 
transforming  is 


Nc=l 

Nu  jo 


(5) 


Since  it  depends  on  both  the  particle  charging  current,  jo,  and 
the  nucleation  rate,  r,  it  provides  a  method  of  showing  how  relative 
changes  in  nucleation  vs.  growth  rates  change  with  charging 
conditions. 


5.  Conclusions 

In  this  study,  we  have  demonstrated  a  number  of  experimental 
techniques  that  shed  insight  on  the  phase  transformation  mecha¬ 
nisms  active  during  LiFePOzi-battery  cycling.  We  describe  a  method 
of  sectioning  electrodes  for  high-spatial  resolution  morphological 
and  chemical  analysis,  which  has  significant  advantages.  We  also 
show  that  the  combination  of  energy-filtered  TEM  and  data  pro¬ 
cessing  routines  that  include  multivariate  statistical  analysis  and 
multiple  least  squares  fitting  can  be  used  to  give  state-of-charge 
maps  for  LiFePC^  electrode  particles.  The  results  are  compared  to 
results  from  the  STXM  and  are  in  good  agreement,  validating  both 
approaches.  The  good  agreement  between  both  datasets  suggests 
that  it  is  possible  to  obtain  the  needed  information  from  the  TEM. 
Since  TEM  is  more  widely  available  than  STXM,  EFTEM  presents  a 
convenient  way  to  probe  nanoscale  chemistry  in  battery  electrodes. 


Finally,  our  results  reaffirm  the  findings  Brunetti  et  al.  [29]  and 
Chueh  et  al.  [4]  that  a  small  fraction  of  mixed  phase  particles  exist. 
Brunetti  et  al.  concluded  that  this  result  provides  support  for  the 
“domino-cascade”  mechanism.  For  the  charging  conditions  of  the 
batteries  examined  in  these  studies,  the  actual  transformation  of  an 
individual  particle  happens  very  fast  relative  to  the  time  it  takes  to 
charge  a  battery.  Therefore,  whether  the  particle  transforms  via  a 
“domino-cascade”  mechanism  or  a  “core-shell”  mechanism  is 
irrelevant  to  optimizing  performance.  The  important  point  is  the 
time  it  takes  for  the  new  phase  to  nucleate.  Because  smaller  par¬ 
ticles  are  more  likely  to  nucleate  the  new  phase  [4],  electrodes  with 
a  smaller  average  LiFePCU  particle  size  may  have  improved  per¬ 
formance.  Alternatively,  some  other  means  of  reducing  the  nucle¬ 
ation  barrier  such  as  introducing  heterogeneous  nucleation  sites 
may  improve  performance. 
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